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We explore the effect of a (non) magnetic impurity on the thermal transport of the spin-1/2 
Heisenberg chain model. This unique system allows to probe Kondo-type phenomena in a prototype 
strongly correlated system. Using numerical diagonalization techniques we study the scaling of the 
frequency dependent thermal conductivity with system size and host-impurity coupling strength as 
well as the dependence on temperature. We focus in particular on the analysis of "cutting-healing" 
of weak links or a magnetic impurity by the host chain via Kondo-like screening as the temperature 
is lowered. 
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I. INTRODUCTION 

A (non) magnetic impurity coupled to a spin-1/2 
Heisenberg chain is a prototype system that exemplifies 
"Kondo"-type effects in a correlated system. Starting 
with the proposal of Kane-Fisher,^ a weak link in a re- 
pulsive (attractive) Luttinger liquid was shown to lead to 
an insulating (transmitting) ground state. The cutting or 
healing of spin chains by a variety of (non) magnetic de- 
fects has also been established^ ^ as well as the effect of a 
magnetic impurity on the ground state of the anisotropic 
easy-plane Heisenberg chain .^^ Generically, a weak link or 
coupling to a magnetic impurity in a Heisenberg antifer- 
romagnetic chain leads to a ground state corresponding 
to two open chains. In the exceptional case of two adja- 
cent links or a ferromagnetic (attractive in the fermionic 
language) easy axis anisotropy a healing of the defect 
is conjectured.^ This screening effect is characterized by 
a Kondo-like temperature and screening length.^ ^ These 
phenomena have so far mostly been studied either as they 
are reflected on ground state properties, e.g. finite size 
gaps, entanglement or, somewhat indirectly, as a temper- 
ature dependent induced staggered susceptibility.l^ 

In this work we use an exceptional physical probe for 
the study of these effects, namely, the thermal trans- 
port in the spin-1/2 Heisenberg chain that is truly 
singular. Although the Heisenberg model describes a 
strongly correlated system, the thermal conductivity is 
purely ballistic as the energy current commutes with the 
Hamiltonian,^ a result that is related to the integrability 
of this model. ^ Thus the only scattering present is due to 
the defect and thus its frequency /temperature/coupling 
strength dependence can be isolated and clearly ana- 
lyzed. In this context it was already found that a 
single potential impurity renders the thermal transport 
incoherent^ with the frequency dependence of the ther- 



mal conductivity well described by a Lorentzian, at least 
for a weak impurity. This is in sharp contrast to the case 
of a non-interacting system where in spite of the impu- 
rity the transport remains coherent described within the 
Landauer formalism by a finite transmission coefficient 
through the impurity. Thus a single static impurity ma- 
terializes the many-body character of scattering states. 

Besides its theoretical interest, the effect of 
(non) magnetic impurities on the thermal transport 
of quasi-one dimensional materials as SrCu02, Sr2Cu03 
and the ladder compound La5Ca9Cu2404i has recently 
become possible to explore experimentally.^^ 

In this work we use numerical diagonalization tech- 
niques - (full) exact diagonalization (ED), the Finite- 
Temperature Lanczos method (FTLM)^^ and the Micro- 
canonical Lanczos method (MCLM)^^ - to study the ther- 
mal transport in the Heisenberg chain model either cou- 
pled to a magnetic impurity or perturbed by single and 
double weak links. These state of the art techniques are 
crucial in the attempt to look for the subtle low tempera- 
ture many-body effects associated with Kondo screening. 

II. MODEL 

We consider the one dimensional anisotropic spin-1/2 
Heisenberg model in the presence of a magnetic impurity 
out of the chain or weak links, 

L-l 
1=0 

hi^i^i = sfsf^^ + 5f 5f+, + A^f^f+i, (1) 

where 5", a = x^y^z are spin-1/2 operators, J^,^+i > 
the in-chain magnetic exchange coupling that we take 
antiferromagnetic, J' the chain-impurity coupling. A, A' 
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anisotropy parameters and S a spin-S magnetic-impurity 
operator {h = 1). In this work we mostly consider a spin- 
1/2 impurity. We assume periodic boundary conditions, 
sl = So, and uniform couplings J^,z+i = J, except in the 
study of weak links (see below). We vary the anisotropy 
parameters A, A', with A = A', in order to look for the 
(healing) cutting of the chain effects mentioned above. 

In our study, based on standard linear response the- 
ory, the frequency uj dependence of the real part of the 
thermal conductivity (regular component) is given by 

o ■ P+OO 

= --X"(^), XH = J / dte^-\[f{t),f]), 
^ ^ Jo 

(2) 

where (3 = 1/T, T is the temperature and ks = 1. 
We determine the energy current from the hydrody- 
namic {q 0) limit of the energy continuity equation 
dHq/dt ^ qf with Hq = e*^^/i/,^+i as, 

L-l 

f = ^ Jl-l,lJl,l+l Si • (s^+i X S^_i) 

JJf 

+ • (S X sl-1 +si X S), (3) 

showing for simplicity the case A = 1 (A 7^ 1 is ob- 
tained by sf Asf in the cross-product terms). When 
J' = and all J/,/+i = J the energy current commutes 
with the Hamiltonian, the transport is purely ballistic 
and the thermal conductivity consists of only a S{uj)-peak 
proportional to the thermal Drude weight. 



III. HIGH TEMPERATURE LIMIT 

Starting from the high temperature {/3 0) limit we 
can obtain a first impression on the behavior of the fre- 
quency dependence of f<i{uj) from the 0th and 2nd mo- 
ments, iin = j dujuj^hz{uj) which are equal to (for the 
isotropic point, A = 1), 

6 2 

/io = const. X ^(J2 + -S^), (const. = tt^) (4) 

/i2 = const. X — — (39J^ - 12JJ' + + 36S^) , 
Ll ^ 

where = {J^ /3)S{S + 1) is the characteristic impu- 
rity spin dependence. One could expect the 2nd mo- 
ment to reflect the width of hz{uj) and thus to be re- 
lated to the inverse scattering time 1/r. We note that 
for this impurity problem an assumption of a Gaussian 
form K,{uj) = /^c^c cxp"^'^'^^ would imply from the L de- 
pendence of /io,2 that tvdc = ^(0) would scale as and 
1/r ^ \l\fL. This is, however, incorrect as is also evi- 
dent from the disagreement with higher moments, n > 2, 
which behave all as \in cx For weak-coupling cases, 
such as a single impurity weakly coupled to the host 



chain, we should therefore rather expect a Lorentzian- 
like frequency dependence with a static cx L and a 
characteristic frequency width 1/r cx 1/L. 
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FIG. 1: Frequency-dependent thermal conductivity in the 
high-T limit scaled as Hi{ujL)/L for (A = 1): (a) weak cou- 
pling J' — 0.5J, (b) strong coupling J' — 2J (curves are 
normalized to unity). 

In Fig. [1] we show the frequency dependence of 
the thermal conductivity, normalized and appropriately 
scaled with system size. Note that in the high-T (/3 0) 
limit the relevant (but still nontrivial) quantity is T'^hz{uj) 
which is implicitely extracted by the normalization. We 
thus present results of the normalized n{ujL)/L for a 
weak, J' = 0.5J and strong, J' = 2J coupling case re- 
spectively. The data up to L = 16 were obtained by 
fuh ED while for L = 18 - 22 the MCLM was used.!^ 
The ^— peaks at the excitation frequencies are binned in 
windows 5uj = 0.01, which also gives the frequency res- 
olution of the spectra. For J' = 0.5 J we find a simple 
Lorentzian form while in the strong coupling case the 
behavior is nonmonotonic with a maximum at a finite 
frequency 0{1/L). In both cases the proposed L scaling 
is indeed realized. 

As for the scaling with impurity spin S suggested 
by the proportionality of the 2nd moment to = 
{r^/3)S{S + 1) we show in Fig. [2] MCLM results for 
hz{ujL)/L for a series of S-values and couplings J' so that 
the effective perturbation strength B^ retains its value. 
We find indeed that at both weak as well as strong cou- 
pling the scaling is well obeyed, giving a wider applica- 
bility to our results. They can be applied to a range of 
impurity spin values becoming directly relevant in the 
interpretation of experiments. 

Now let us address the generic L ^ 00 behavior. We 
can discuss it by considering the memory function N{uj) 
representation defined via the general complex function 

Xo 

^(^) = ^ ATf \ ^ Xo=x(^^0). (5) 

UJ -\- I\ yUJ ) 

where the real hz{uj) = R' {uj) and N"{uj) ~ 1/r plays the 
role of the (frequency dependent) thermal-current relax- 
ation rate. The lowest moments jin can be evaluated (in 
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b) J' S(S+1)«1.70 




5 10 15 20 

ooL 



FIG. 2: Frequency dependence of the normalized thermal 
conductivity k{ujL)/L in the high-T limit for a variety of 
impurity spin values S = 1/2,1,3/2,2 and for: (a) J' / J = 
0.5,0.3,0.22,0.18 corresponding to the weak coupling = 
'^/3)S{S + 1) ^ 0.06, (b) f/J = 1.5, 0.92, 0.67, 0.53 corre- 
sponding to the stronger coupling 0.57. 




FIG. 3: Memory function (uj) for a strong coupling J' = 2J 
and for various lattice sizes L — 12 — 24, using both ED and 
FTLM. Inset: the scaled function N"{ijoL) is shown at low 
frequencies. 



principle) exactly in the high-T limitP^on a finite size lat- 
tice of L sites. Involving only local quantities, at least for 
< n < L/2, they should behave as fin = fin/L whereby 
jin is size independent for n < L/2. It is plausible that 
also higher moments, n > T/2, behave as /i^ oc 1/L. If 
/in for n > L/2 would be also size independent, then this 
would imply the scaling N{u) = -^7V(cj), with a universal 
(size independent) N{uj). Consequently 



{ujL) + N{uj) ' 



with the real part i<i{(jj) for L oo and u 
the Lorentzian scaling relation. 



(6) 



obeying 



k{loL) ^ l3xoN"iuj 
L 



0) 



(wL)2 + Ar"(w ^ 0)2 



(7) 



provided that N^'{uj ^ 0) is finite. This is, however, 
clearly not what we observe in Fig. [T] where from the 
non-Lorentzian shape we must conclude that the memory 
function also scales as N{ujL) and thus. 



hz{ujL) 
L 



(8) 



This is not in contradiction with the moments argument, 
since the higher moments, n > T/2, determine the low 
frequency behavior. So we can argue that at high fre- 
quencies N{uj) scales as u while at low frequencies as ooL. 
This scenario is indeed verified in Fig. [sjat the low/high 
frequency regimes, where N{uj) is extracted from the 
hz{uj) data. The FTLM method is used for lattice sizes 
L > 16 with Ml = 500 Lanczos steps and smoothed 
with an additional frequency broadening Suj = 0.03. On 
the other hand, we can also explain the observed general 
K:{ujL)/L scaling with the similarity to a noninteracting 
system - with an impurity. In the latter case, the char- 
acteristic scaling Luj is signature of "free" oscillations in 
the system. 
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FIG. 4: Impurity coupling J' dependence of scaled 
N"[Luo)/J''^ and the comparison with the perturbative result. 
Results are obtained for A = 1 and L = 16 via ED. 

To study the crossover from weak to strong coupling 
regime we show in Fig. |4]the evolution of the relaxation- 
rate function N"{uoL) with impurity coupling J' along 
with a perturbative evaluation Nq{ujL) using the eiges- 
tates of the Hamiltonian without the impurity. It is in- 
teresting that the memory function shows an increasingly 
pronounced structure with minima at approximately the 
same frequencies, multiples of 27r/L independently of J' 
and which are not present in the perturbative calculation. 
In particular the characteristic frequency of the minima 
decreases as the anisotropy parameter A decreases and 
thus it apparently related to the velocity of elementary 
excitations (spinous) in the system. We can conjecture 
that this peak structure is due to a resonant mode, cre- 
ated by multiple forward/backward scattering on the im- 
purity, characteristic for the noninteracting system. It is 
remarkable that this happens even in this high tempera- 
ture limit. This effect has already been seen in integrable 
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systems where a perturbation seems to affect the totahty 
of the energy spectrum. Now the picture is clear, N^\uj) 
increases as J'^, scales as ujL at low frequencies and at 
the same time develops a structure that dominates the 
behavior of KiujL) turning the Lorentzian weak-coupling 
shape to a nontrivial one at strong coupling. 



IV. WEAK LINKS - FINITE T 

Next we examine the behavior of the thermal conduc- 
tivity n{uo) as we lower the temperature, starting with 
the influence of static weak exchange links. 

Kane-Fishei^ for a Luttinger liquid and Eggert and 
AflleclP (EA) for the isotropic spin- 1/2 Heisenberg chain, 
proposed that a weak link leads to an open chain (cut- 
ting) in the low energy limit. In contrast, a defect of 
two adjacent weak links is "healed" leading to a uniform 
chain at T = 0. 
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FIG. 5: Frequency dependence of: (a) the normalized ther- 
mal conductivity k{ujL)/L^ (b) the extracted memory func- 
tion N" (ujL)^ for a chain of L = 22 sites with one weak link 
J = 0.7 J and various Tj J — 0.3 — 2.0 . (c) Temperature 
dependence of Hidc{T)/L. 

To analyze this effect we consider a chain with only 
one weak link, that is one altered bond with coupling 
e.g. Jo,i = J in an otherwise uniform chain {J' = 0, 
there is no spin impurity). The characteristic Kane- 
Fisher temperature is given in the weak coupling limit 
by Tkf ^ {J — JY / J ' In Fig. [s^i we show the corre- 
sponding f<i{ujL)/L for J = 0.7 J and a series of temper- 
atures. The data are obtained using the FTLM method 
for a chain of L = 22 spins, by Ml = 2000 Lanczos 
steps and smoothed by an additional frequency broaden- 
ing 5uj = 0.007. From Fig. [S^i we notice that t<i{ujL)/L 
develops a strongly nonmonotonic frequency dependence 
by lowering the temperature, with a maximum at a flnite 
frequency that suggests a flow to the strong coupling limit 
similar to the one discussed before by increasing J' . In 



Fig.jsjD, the extracted N"{ujL) for various T is presented, 
with the development of a characteristic structure that 
explains the nonmonotonic behavior of hi{uj). The in- 
creasing value of A^''(0) ~ 1/r with decreasing tempera- 
ture indeed corresponds to the effect of "cutting" of the 
chain. 

Nonmonotonic is also the frequency dependence of 
hi{ujL)/L for the case of two adjacent equal weaker links, 
Jl-1,0 = Jo,i = J = 0.7J, as shown in Fig. [i^. How- 
ever, in this case we observe in Fig. [gJd the opposite 
behavior of N^\uj). Namely "healing" of the double de- 
fect deduced by the decreasing A^''(0) as the temperature 
is lowered in agreement with theoretical prediction.^ We 
should note that both cutting/healing are low frequency 
effects at frequencies ujL 0(1). 
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FIG. 6: Frequency dependence of: (a) the normalized ther- 
mal conductivity k{(jjL)/L, (b) the extracted memory func- 
tion N"(uoL) for a chain of L = 22 sites with two adjacent 
weak links J = 0.7J and various T/J = 0.3 - 2.0. (c) Tem- 
perature dependence of Hidc{T)/L. 
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FIG. 7: Temperature dependence of iV''(0) for J = 0.5,0.7J 
showing cutting/healing behavior for one and two weak links. 

To summarize the observed behavior we show in Fig. 
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j?! the T-dependence of the relaxation rate N"{0) for two 
different couphngs J/J = 0.5,0.7, for one and two weak 
hnks, respectively. The presented results confirm the ex- 
istence of the cutting behavior at low T for a single link, 
as well as the healing by lowering T for two adjacent and 
equal links. As expected, both effects appear only at low 
T/ J < 1 while the dependence of the characteristic Tkf 
on J/ J is less pronounced. 



V. SPIN COUPLED TO THE CHAIN - FINITE T 

Finally we can study the effect of lowering the tem- 
perature on the scattering by a magnetic impurity. Ac- 
cording to EA it leads to cutting the chain at T = ir- 
respective of the sign of J^ This proposal was extended 
by Furusaki and HikiharsP to the anisotropic spin chain 

— 1 < A < 1 where they furthermore proposed that for 

— 1 < A < (attractive case in the fermionic language) 
there is "healing" of the impurity, in analogy to the case 
of two adjacent weak links. 

In the Kondo problem the characteristic temperature 
in the weak coupling limit is given by Tk ~ v exp(— c/ J') 
with c being a constant, v the velocity of spin excitations 
and J' the Kondo coupling. In the case of a spin-1/2 
chain it was showiP^ that th e exponential depe ndence 
is replaced by Tk - exp(-7rv'l/J' - {S' + 1/2)2) g^^d a 
next-nearest neighbor coupling J2 ^ 0.2412 is needed to 
recover the traditional Kondo case. We should note that 
in the model studied the impurity spin is attached only 
at the end of the chain - in contrast to our model - but 
plausibly the behavior is qualitatively similar. To get a 
qualitative idea of orders of magnitude for our proble nP 
for J' = 0.3J, Tk - 0.014, ^k ^ 40, for J' = 0.6J, 
Tk - 0.388, ^k ^ ^ and J' = J, ^k = 0.65. As in 
our study we are limited to T > 0.4 in order to see a 
"Kondo" crossover we must consider a coupling J' > 0.5 J 
and thus we are in the relatively strong coupling regime, 
with typical screening length of the order 1- 

In Fig. 11 we show n{(jjL)/L for a chain of L = 22 
sites at strong coupling J' = 2J and two representative 
cases A = ±0.5 as we lower the temperature. Indeed 
we find at low frequencies the gradual development of 
the corresponding "cutting/healing" behavior which we 
exemplify in the inset by A^''(0) as a function of temper- 
ature both for A = ±0.5 and the most typical isotropic 
case A = ±1.0. It is remarkable that the tendency to 
increase-decrease the scattering time is already evident 
from high T, presumably due to the local character of 
the effect because of the strong J' coupling. We note in 
passing that the ujL scaling is found not just at high T 
but rather at all T (not shown). 

Next in Fie;. [9] we show N"{{)) as a function of T for 
a series of increasing J' couplings. The "cutting" effect 
for the repulsive case A = ±0.5 is present for all values 
of J' with no easily distiguishable "Kondo" temperature. 
We are always dealing with screening lengths well less 
than the system size where presumably no subtle many- 
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FIG. 8: Frequency dependent normalized thermal conductiv- 
ity k{ujL)/L for strong coupling J' = 2 J, A = ±0.5 and three 
T/J = 50, 2, 0.4. Inset: T-dependence of N" {0) for A = ±0.5 
and A +1. 



body effects come into play. On the other hand, in the 
attractive case A = —0.5, we do not observe "healing" 
for the weakest coupling J' = ±0.5 where the screening 
length is expected to be several lattice sites. 
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FIG. 9: A^^'(O) vs. T for the repulsive (attractive) case A 
±0.5(-0.5) for different J' / J = 0.5, 1.0, 1.5. 



Finally, in Fig. [To] we summarize the T-dependence 
of ndc/L for a variety of coupling strenghts J' / J and 
A = ±0.5. The experimentally most interesting case 
A = ±1 corresponding to isotropic antiferromagnetic as 
well as ferromagnetic impurity coupling is shown in Fig. 
For A > we observe in Fig. 
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10 i and Fig. 11 



continuous decrease of the Kdc with increasing J' . This 
can be explained with the formation of a local singlet, 
at least for T < J' which blocks the current through the 
impurity region. On the other hand, the A < case in 
Fig. P!q) 3 reveals a saturation of i^dc with J', at least for 
intermediate large J' . However, for severe perturbations 
{J' ^ J) the impurity cannot be healed by the chain 
leading inevitably to a further decrease of the Kdc- 
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FIG. 10: Temperature dependence of Kdc/L for a variety of 
impurity couplings J' and for: (a) repulsive A — +0.5, (b) 
attractive A = —0.5. 
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FIG. 11: Temperature dependence of tvdc/L for a variety of 
impurity couplings A = 1 and for: anti- ferromagnetic 
couplings (top), ferromagnetic couplings (bottom). 



VI. CONCLUSIONS 
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FIG. 12: Frequency dependence of the thermal conductivity 
k{(jjL)/L in the high temperature limit for various values of 
the coupling f/J = 0.8 - 4.0 and A = 1.0. 



{b) Furthermore, universal oscillations in the dynamical 
relaxation rate N^\uj) become visible, from the weak cou- 
pling regime already, with the period uj (x 1/L being a 
remnant of the impurity multiple-scattering phenomena 
in a noninteracting system. 

(c) Our results confirm the existence of the Kondo-type 
effects of impurities on lowering the temperature. In the 
case of weak links and for the isotropic Heisenberg model 
cutting and healing effects are observed at lower T for a 
single weak link and a pair of identical weaker links, re- 
spectively, in accordance with theoretical predictions.!^ 
In the case of a spin coupled to the chain the cut- 
ting/healing effects at low T depend on the sign of the 
anisotropy A. For ferromagnetic anisotropy (A < 0), 
the chain screens the impurity and the system enters the 
weak coupling regime as the temperature is decreased. 
The opposite behavior is obtained for antiferromagnetic 
anisotropy (A > 0), where the system flows to the strong 
coupling limit at lower temperatures. 
{d) Obtained data can be used to model the behavior 
observed in experiments on materials with spin chains 
doped with magnetic and nonmagnetic impurities.^ 



In conclusion, by analysing the unique behavior of the 
thermal conductivity of the spin- 1/2 Heisenberg model 
several effects of the local static and dynamical impuri- 
ties have been established: 

(a) A single local impurity, either static as the local field^ 
and weak link, or dynamical as the spin coupled to the 
chain turn the dissipationless thermal conductivity into 
an incoherent one. Numerical results for the dynamical 
conductivity, best studied at high-T, reveal that a single 
impurity in a system of L sites shows a universal scaling 
form f<i{ujL)/L at least in the low-o; regime. For weak 
perturbation, as weakly coupled spins outside the chain, 
the scaling form is of the simple Lorentzian type. On the 
contrary large local perturbation can lead to a nontrivial 
form with the maximum response at > 0. 
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Appendix A: Open chain 

Throughout the article the term "cutting" is used to 
describe the behavior of the system in the strong coupling 
limit. In order to justify the term "cutting", we present 



in Fig. 12 results for the thermal conductivity of a chain 
of L = 16 sites obtained by ED in the high temperature 
limit for various couplings J' and the thermal conduc- 
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tivity of a uniform chain with open boundary conditions 
as weh. Fig. 12 ihustrates the flow of the system from 
a Drude hke behavior (weak couphng) to a chain with 
open boundary conditions (strong couphng), which was 
aheady proposed for a single non-magnetic impurity (a 
local field) from the level statistics analysis.^ We choose 
to present the jagged results, i.e. without implement- 
ing any smoothing procedure, in order not to wash out 



the development of the narrow peaks corresponding to 
the excitations of the open chain. For the strong cou- 
pling cases there is some rather significant structure at 
frequencies ~ J' which correspond to local excitations 
of the impurity. However, these excitations are irrele- 
vant for the effect of the impurity on the chain which is 
studied here. 
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